Spin effects on gravitational waves from inspiraling compact binaries at second 

post-Newtonian order 
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We calculate the gravitational waveform for spinning, precessing compact binary inspirals through 
second post-Newtonian order in the amplitude. When spins are coUinear with the orbital an- 
gular momentum and the orbits are quasi-circular, we further provide explicit expressions for 
the gravitational-wave polarizations and the decomposition into spin-weighted spherical-harmonic 
modes. Knowledge of the second post-Newtonian spin terms in the waveform could be used to 
improve the physical content of analytical templates for data analysis of compact binary inspirals 
and for more accurate comparisons with numerical-relativity simulations. 
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I. INTRODUCTION 

Coalescing compact binary systems are a key source 
of gravitational radiation for ground-based gravitational- 
wave detectors such as the advanced Laser Interferom- 
eter Gravitational wave Observatory (LIGO) the 
advanced Virgo [1], the GEO-HF [1], the Large Cryo- 
genic Gravitational Telescope (LCGT) (or KAGRA) jll, 
coming into operation within the next few years, and 
future space-based detectors 0, Q. For this class of 
gravitational-wave sources, the signal detection and in- 
terpretation will be based on the method of matched fil- 
tering 0, Q : where the noisy detector output is cross- 
correlated with a bank of theoretical templates. The ac- 
curacy requirement on the templates is that they remain 
as much as possible phase coherent with the signal over 
the hundreds to thousands of cycles of inspiral that are 
within the detector's sensitive bandwidth. 

Constructing such accurate templates has motivated a 
significant research effort during the last thirty years. In 
the regime where the separation between the two bodies 
is large, gravitational waveforms can be computed using 
the post-Newtonian (FN) approximation method 
In the post-Newtonian scheme, the results arc written 
as an asymptotic expansion in powers of va/c, with va 
being the magnitude of the orbital coordinate velocity 
Va of body ^ at a given time. This approximation is 
physically relevant for va/c ^ 1, i.e. in the so-called 
inspiraling regime where the radiation reaction forces, of 
order (va/c)^ are negligible over an orbital period and 
act adiabatically on a quasi-conservative system. In the 
domain of validity of the post-Newtonian scheme, the 
separation r ^ (GniA/ c^) ^ ("^/c)^, with m ~ mi + m2 
and V = \v\ = \vi — V2\, remains small with respect to 
the radii of both compact objects ~ GmAj (? or, in other 
words, the bodies can be regarded effectively as point 
particles. 

Post-Newtonian waveforms cease to be reliable near 



the end of the inspiral and the coalescence phase, where 
numerical-relativity simulations should be used to pre- 
dict the gravitational- wave signal By combin- 
ing the information from post-Newtonian predictions and 
the numerical-relativity simulations it is possible to ac- 
curately and analytically describe the gravitational- wave 
signal during the entire inspiral, plunge, merger and ring- 
down process [l5l - l23j . 

For nonspinning binaries, the post-Newtonian expan- 
sion has been iterated to 3.5PN order beyond the leading 
Newtonian order in the gravitational-wave phasing [2J- 
[26j . The gravitational- wave amplitude has been com- 
puted through 3PN order 



mode through 3.5PN order 



27-|30| and the quadrupole 
3l|. However, black hole bi- 
naries could potentially have large spins j32| which may 
be misaligned with the orbital angular momentum, in 
which case the precession effects add si gnif icant complex- 
ity to the emitted gravitational waves [33|. Ignoring the 
effects of black hole spins could lead to a reduction in 
the signal-to-noise ratio and decrease the detection ef- 
ficiency [33, HH although this should be overcome with 
phenomenological and physical models Hi], [36l - l43j . To 
maximize the payoffs for astrophysics will require extract- 
ing the source parameters from the gravitational-wave 
signal using template models computed from the most 
accurate physical prediction available [43 - l47| . Spin ef- 
fects in the waveform are currently known through much 
lower post-Newtonian order than for nonspinning bina- 
ries. More specifically, s pin effects are known through 
2.5PN order in the phase |48l- [50| . 1.5PN order in the po- 
larizations for spin-orbit effects j5ll[52| . 2PN order for the 
spini-spin2 effects [5l|,|5l] and partially 3PN order in the 
polarizations for the tail- induced spin-orbit effects (53 |. 

In this paper, we compute all spin effects in the 
gravitational- wave strain tensor through 2PN order. This 
requires knowledge of the infiuence of the spins on the 
system's orbital dynamics as well as on the radiative mul- 
tipole moments. At this FN order, nonlinear spin effects 
due to the spin-induced quadrupole moments of the com- 
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pact objects first appear. Using results from Ref. [55l - l58| . 
we derive the stress-energy tensor with self-spin terms 
and compute the self-induced quadrupole terms in the 
equations of motion and in the source multipole moments 
at 2PN order. Our results are in agreement with previous 
calculations j59l - [6^ . 

The two main inputs entering our calculation of the 
gravitational-wave strain tensor through 2PN order are 
(i) the results of Refs. [13, HH HI] for the influence of 
the spins on the system's orbital dynamics, which have 
also been derived by effective field theory and canonical 
methods (sgI. l63l - [67| . and (ii) the spin effects in the sys- 
tem's radiative multipole moments [H^l- Recently, the 
necessary knowledge to compute the waveform at 2.5PN 
order was obtained using the effective field theory ap- 
proach [i^ll^. Here we use (i) and (ii) in the multipolar 
wave generation formalism |68l47C| to obtain the wave- 
form for spinning, precessing binaries through 2PN order. 
To compute the gravitational polarizations from this re- 
sult, one must specify an appropriate source frame and 
project the strain tensor onto a polarization triad. For 
precessing systerns, there are several frames that could be 
employed d, [H, [5l|, [Tll - lTSj . For non- precessing binaries 
with the spins collinear to the orbital angular momen- 
tum, the most natural frame is the one used for nonspin- 
ning binaries. Therefore, instead of choosing one frame, 
for simplicity, we specialize to the non-precessing case 
and quasi-circular orbits and provide the explicit expres- 
sions for the gravitational polarizations. Lengthy cal- 
culations are performed with the help of the scientific 
software Mathematica®, supplemented by the package 
xTensor [tI] dedicated to tensor calculus. Our generic, 
precessing result is available in Mathematica format upon 
request and can be used to compute the polarizations for 
specific choices of frame. We notice that the 2PN terms 
in the polarizations, for circular orbits, linear in the spins 
were also computed in Ref. [t^]- However, these results 
contain errors in the multipole moments, which were cor- 
rected in Ref. jsof . 

For future work at the interface of analytical and nu- 
merical relativity, we also explicitly compute the decom- 
position of the strain tensor into spin- weighted spherical- 
harmonic modes for non-precessing spinning binaries on 
circular orbits. The test-particle limit of these results can 
also be directly compared with the black-hole perturba- 
tion calculations of Refs. [zl, [t^, and we verify that the 
relevant terms agree. 

The organization of the paper is as follows. In Sec. [Hi 
we review the Lagrangian for compact objects with self- 
induced spin effects [ssi - fSTl |6]| . compute the stress- 
energy tensor and derive the self-induced spin couplings 
in the two-body acceleration and source multipole mo- 
ments ^59t - i62| . In Sec. IIIII we summarize the necessary 
information about spin effects in the equations of motion 
and the wave generation necessary for our calculation. 
In Sec. IIVBI we calculate the spin-orbit effects at 2PN 
order in the strain tensor for generic precessing binaries. 
In Sec. IIV CI we complete the knowledge of 2PN spin- 



spin terms by including the spin self-induced quadrupole 
terms in addition to the spini-spin2 terms obtained in 
Ref. |5l|. In Sec. IIV El we specialize to quasi-circular or- 
bits and explicitly give the polarization tensors for non- 
precessing systems. Then, in Sec. IIVFI we decompose 
the polarizations into spin-weighted spherical-harmonic 
modes. Finally, Sec. IVl summarizes our main findings. 

We use lowercase Latin letters a, i, j, ... for in- 
dices of spatial tensors. Spatial indices are contracted 
with the Euclidean metric, with up or down place- 
ment of the indices having no meaning and repeated 
indices summed over. We use angular brackets to de- 
note the symmetric, trace-free (STF) projection of ten- 



sors, e.g., T( 



STF[r,j 



\&ijTkk, where 



the round parentheses indicate the symmetrization op- 
eration. Square parenthesis indicate antisymmctrized in- 
dices, e.g., Tjjj] = \{Tij — Tji). The letter L = 
signifies a multi-index composed of i STF indices. The 
transversc-traceless (TT) projection operator is denoted 
^ijlh = 'Pa{^VJ)b-^'P^JVab, whcrc T'y = S^j-N,Nj is the 
projector orthogonal to the unit direction N = X /R of 
a radiative coordinate system X'^ = {cT,X), where the 
boldface denotes a spatial three-vector. As usual, 5^1, 
represents the space-time metric and g its determinant. 
The quantity Sijk is the antisymmetric Levi-Civita sym- 
bol, with £123 = 1, and e^^po- stands for the Lcvi-Civita 
4-volume form, with eoi23 = +^/—g- Henceforth, we shall 
indicate the spini-spin2 terms with 5152, the spinf , spin^ 
terms with 5^ and the total spin-spin terms with SS. 
Throughout the paper, we retain only the terms relevant 
to our calculations and omit all other terms, which are ei- 
ther already known or at a higher post-Newtonian order 
than required for our purposes. 



II. MODELING SPINNING COMPACT 
OBJECTS WITH SELF-INDUCED 
QUADRUPOLES 

In this section we review the construction of 
a Lagrangian for compact objects with self-induced 
quadrupole spin effects jssl - ISTl [61I Isoj . compute the 
stress-energy tensor and derive the self-induced spin cou- 
plings in the two-body acceleration and source multipole 
moments. Our findings are in agreement with previous 
results ii-il. 



A. Lagrangian for compact objects with 
self-induced spin effects 

A Lagrangian for a system of spinning compact objects 
with non-dynamical^ self-induced quadrupole moments 



^ We shall not include kinetic terms in the Lagrangian for the 
quadrupole moment that can describe resonance effects in neu- 
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can be obtained by augmenting the Lagrangian for point 
particles with describing the quadrupole-curvature 
couphng for each body A. Since the action for body A 
must admit a covariant representation, the corresponding 
Lagrangian should be a function of the 4-velocity 



M^, the metric p^i,, 



the Riemann tensor R and its 



covariant derivatives, evaluated at the worldline point 
j/^, and the spin variables entering via the antisymmetric 
spin tensor S'^ . 

The spin tensor S^' contains six degrees of freedom. 
It is well known that in order to reduce them to the three 
physical degrees of freedom a Spin Supplementary Con- 
dition (SSC) should be imposed [8l|. This is equivalent 
to performing a shift of the worldline y^. In this paper 
we specialize to the SSC Sj^p^ = which is equivalent 
to S'^u^ = since « ttiacu^ through 2.5PN order. 
To ensure the preservation of the SSC under the evolu- 
tion, we follow Ref. [53| and introduce the spin tensor 

The spin tensor S'T' automat- 



s7- 



'A " 



ically satisfies the algebraic identity S' 



0. which 



provides three constraints that can be used to reduce the 
spin degrees of freedom from six to three. 

From the above discussion and Refs. [H, [13, we as- 
sume that the Lagrangian of particle A is of the form 



^A 



LAtivXpS^ S/ , where LAtivXp is a polynomial in 
the Riemann tensor and its derivatives, as well as the 
4-velocity u^. As noticed in Ref. [s^, any term propor- 
tional to ^ ...Rafj evaluated at point can be recast 
into a redefinition of the gravitational field. As a re- 
sult, the Riemann tensor may be replaced for each of 
its occurrences by the Weyl tensor C^p^^, which can be 
decomposed into a combination of the gravito-electric 
and gravito-magnetic type STF tidal quadrupole mo- 



^3 T}A* 



pa 
af}- 



if R'^ 
2 '^pvpa-^^ 



X 
ppv 



ments G ^ 
HpAVa) ^ 

More generally, the multiple space derivatives of C 
at point j/^ may be expressed in terms of some STF tidal 
multipole moments and i?^^ . . .^^, of parity 1 and 

— 1 respectively. However, those higher-order moments 
will play no role in this paper. 

Taking into account that the contraction of the veloc- 
ity vector with both G^^, and S'^ vanishes, that the 
spin and tidal multipole tensors are traceless, and that 
the Lagrangian must obey parity and time-reversal sym- 
metries we obtain I55l-l57l [831 



(2.1) 



Here, we have also assumed that the rotating body is ax- 
ially symmetric and we have replaced Sj^ with Sj^ since 
the difference between these spin variables contributes to 

the equations of motion at 0(5^), where Sa — \J\SaS^\ 

with = £pa,^p.S'7pA/{'2niAc). 



For a neutron star the numerical constant ka in 
Eg (|2.ip depends on the equation of state of the fluid 
[84 | . For an isolated black hole = 1 [H, [fiO] , but for 
a black hole in a compact binary ka can deviate from 1. 
However, these deviations occur at PN orders that are 
much higher than the ones considered here. We notice 
that the leading contribution ka — ^ can be obtained by 
computing the acceleration of body A from Eq. (|2.ip in 
a compact binary for ttia <C m and match it with the 
acceleration of a test particle in the gravitational field of 
a Kerr black hole of mass m (831 . 



B. Effective stress-energy tensor with self-induced 
quadrupoles 

The piece of the stress-energy tensor encoding the self- 
induced quadrupole dynamics of body A reads by defini- 
tion 



rppV 

quad, A 



dTALf[y^{rA),Sf{TA)], 



(2.2) 

where is the Lagrangian (|2.ip . In order to deter- 
mine the action of the operator S/Sg^i,, which stands 
for the usual "functional derivative" with respect to the 
field 5^1, , we need to adopt a specific model for the spin. 
The rotational state of the extended object A is usually 
represented by a tetrad of orthonormal vectors e^— (r^) 
with a e {0,1,2,3} along the worldline with affine 
parameter ta- The corresponding angular rotation ten- 



r^'^Pe%^De\JdTA. We now 



sor is then defined as — // <^Aa" " Ap' 
make the reasonable physical hypothesis that the rota- 
tion of the axially symmetric object takes place about 
the symmetry axis. The moment of inertia I a along that 
direction is a 2PN-order quantity ^ G^m\/c* for com- 
pactness parameters of order 1, whereas Q'^ ^ Va/ Ra, 
Ra being the radius of body A and Va its typical inter- 
nal velocity, is roughly equal to / {Gitia)- In the weak 
field limit where G goes formally to zero, the spin must 



satisfy the relation S"^ = Ia^a ' special relativ- 

ity (85|. In the presence of a non- negligible gravitational 
field, this relation is expected to be modified by non- 
minimal coupled terms proportional to positive powers 



of i?^„„« times positive powers of I a and S'^ 



Ia 



(2.3) 



Here we use a hat to distinguish the generic spin variable 
from the one related to our specific spin model. The 
corrections Ia x 0{Sa/c^) are not relevant for the two- 
body dynamics in this paper because corresponds to the 
4.5PN order when taking into account the factor 0(l/c) 
contained in the spin variable. 

Using the definition (|2.3p for the spin variables, we 
compute in a covariant manner the variation of the action 



tron stars. 



= / dTAL'X {ta) 
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when the metric varies by 5g^y{x), and find tlie following 
quadrupolar piece of the stress-energy tensor 



quad , A 



^2n(A' „'')cA C<^P„.T < r<(f^ q^i c 



AX'-'Ap'-'A 



2VaVp 



(2.5) 



where we have indicated with the Dirac-type scalar 
density = / dA^ S^x'' - y'XiXA)) I V^g(x^ and, 

in the last term, we have adopted the convention that 
symmetrization of indices applies after antisymmetriza- 
tion. As derived in Ref. [80|, the most general form of 
the effective stress-energy tensor is 



Consistently with the approximation already made in the 
spin model (|2.3p . we have neglected here the difference 



(2.4) of order Ia x 0{Sa/c^) between the spins S^j^ and 5^ 



in the above formula. The net result is that Eq. (|2.9ap 
matches Eq. (|2.6p for I = 0, L 2 as expected. Moreover, 
Eq. ((2:9b| agrees with Ref. 

Lastly, the conservation of the stress-energy ten- 
sor (|2.9aP is equivalent to the equation of motion for the 
particle worldline, supplemented by the spin precession 
equation [sl]. They read 



dTA 



ApvX'-' A ' 

(2.10a) 



(2.10b) 



Those equations are in full agreement with the equations 
of evolution derived from the Dixon formalism truncated 
at the quadrupolar order [sgj ]. 



TZi.a{^') 

+ 00 



^Va,Va,...V; 



iTA)n*A{xn , (2.6) 



e=o 



where ta is the proper time of the A*'^ worldline at event 
y'^ with = x° and the coefBcients t'^^^^^^'"'^' {ta), 
are the "skeleton" multipole moments. The latter are 
not arbitrary but satisfy algebraic constraints imposed by 
the equation of conservation VjyTj^^j = 0. Let us check 
that we can indeed recast the total stress-energy tensor, 
including the monopolar, dipolar and quadrupolar pieces, 
in the form ()2.6|) . If we add T^^^^^to the monopolar and 
dipolar contributions jiol I80l. l86l488j 



moii-t-dipole, A 



^ \n\p^^u2c + Vx(^n*ACU'^S'2^ 



(2.7) 



and redefine the spin variable entering the quadrupolar 
piece, as 



ntJ,v ntJ,v 2k^ q\[pi^v\ 

^A - '^A ,^2-'^'^ A '-'AX ' 

we obtain the total stress-energy tensor in the form 



(2.8) 



* ^2 jKiJ-^)p 

A 



(2.9a) 



where the 4-rank tensor J^^'^ takes the following expres- 
sion in our effective description: 



jXppv _ "il^A qa[X p\ qA[^I.„y\ 

2 A "A'-'o- "'A 



mAC 



(2.9b) 



C. Self-induced quadrupole terms in the 2PN 
binary dynamics and source multipole moments 

Once the stress-energy tensor has been derived, the 
post-Newtonian equations of motion and the source 
multipole moments parametrizing the linearized gravi- 
tational field outside the system can be computed by 
means of the usual standard techniques [l^]- At 2PN 
order, the accelerations including the self-spin interac- 
tions were obtained in Refs. [s^, |60|, but the self- induced 
quadrupole effects in the source multipole moments were 
never explicitly included in the standard version of the 
post-Newtonian scheme, although recently they were cal- 
culated at 3PN order using effective-field-theory tech- 
niques [s^l- Here we can use the results of the previous 
section, which constitutes a natural extension of the stan- 
dard post-Newtonian approximation for spinning com- 
pact bodies [i^, and explicitly derive the self- induced 
quadrupole couplings in the 2PN dynamics and source 
multipole moments. 

Henceforth, we define the spin vectors S\ by the re- 
lation Sf/c — g:^Sj^, where Sf is the 3-form induced 
on the hypersurface t = const, by S^. Note that it is 
S\ I c that has the dimension of a spin, while S\ has been 
rescaled in order to have a non-zero Newtonian limit for 
compact objects. 

In the post-Newtonian formalism for point particles in 
the harmonic gauge, it is convenient to represent effec- 
tively the source by the mass density a = (r°°-|-r")/c^, 
the current density a; = T'^'^/c, and the stress density 
(Jij = T'-' . They are essentially the components of the 
stress-energy tensor rescaled so as not to vanish in the 
formal limit c — for weakly stressed, standard matter. 
At 2PN order, the second term in the right-hand side of 
Eq. (|2.9ap does not contribute. From the last term, we 
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obtain the following self-spin contributions: 



+ 1^2 + 0{^ ] , (2.11a) 



' S A 



4- = ^^) 



(2.11b) 
(2.11c) 



where 10 2 represents the counterpart of the preceding 
term with particle 1 and 2 exchanged, and 5i = S'^{x — 

yi). 

At 2PN order, the spin^ part of the equations of mo- 
tion p.lOap for, say, the first particle, reduce to 



cLti 



non-Sl terms - ^dkRloj^S{^ S'^ +0(^) 



(2.12) 

The only occurrence of self-spin interactions at 2PN order 
on the left-hand side of the above equation comes from 
the gradient of the time component of the metric, goo = 
-1 + 2V / + 0{l / c^), where the Newton-likc potential V 
satisfies UV = —AnGa. Although V coincides with the 
Newtonian potential U in the leading approximation, it 
contains higher order corrections, including quadratic-in- 
spin terms coming from the mass density (|2.11ap . which 
are smaller than J7 by a factor 0{l/c^). They read: 



v.: 



2ttGki 



6iyA-i[(5iS'fS'i^^'] + 1 2 



= ^^d,,-Sfs'^+l^2 + o(^), (2.13) 
2mic"^ ri V / 

with di = d/dx' and ri = \x — yi\, the symbol hold- 
ing for the retarded integral operator. Other potentials 
appear at the IPN approximation or beyond, but their 
sources cannot contain a self-induced quadrupole below 
0{\/c^)] thus they are negligible here. The self-induced 
spin part of the acceleration Oi of the first particle is 
therefore given by 



2m{ 

(2.14) 

Replacement of the Christoffel symbols T^^^, and the Rie- 
mann tensor by the leading order values 



with U = Gmi/ri + Gm2/r2 + 0{l/c^ 
explicit result (posing du = d/dy\): 



(2.15) 

yields the more 



1 r '*2 



7712 



Qj ok 



m2Ki 



m^ 



(2.16) 



which agrees with Refs. [5^, [601 '^^ the center-of-mass 
frame, for S^Jc = e.jkS^'' + 0{l/c^). 

Self-induced quadrupolar deformations of the bodies 
also produce 2PN-order terms in the source multipole 
moments 1^ and Jl- Those are defined as volume in- 
tegrals whose integrands are certain polynomials in the 
densities cr, Gi and well as some gravitational po- 

tentials, such as V, that parametrize the metric. Now, 
since those potentials are multiplied by prefactors of or- 
der 0{l/c^) and cannot contain themselves spin^ interac- 
tions below the 2PN order, monomials involving one po- 
tential or more may be ignored for the calculation. The 
remaining sources are linear in the a variables. With the 
help of the general formula (5.15) of Ref. [9l|, it is then 
immediate to get the self-spin contribution to 1^: 



If 



(2.17) 



Inserting the expression p.llap for (Tg2 and performing 
a straightforward integration, we arrive at 



,,„ij,a.,fai'..„;'')sfs« + i"2 + o(£i). 

(2.18) 

We can show similarly that Jl is of order 0{S'\/c'^). 
As a result, at the accuracy level required for the 2PN 
waveform, the only terms quadratic in one of the spins 
that originate from the source moments come from the 
quadrupole £ = 2, for which we have 



102 + 0(1 



(2.19) 



whereas similar terms in {lL)e>3 or {Jl)i>2 lie beyond 
our approximation. The above correction to the mass 
quadrupole agrees with that of Porto et al. jo^l trun- 
cated at 2PN order. It is formally of order ©(l/c"*) but, 
because Sa = ©(l/c^), it is cast to the 3PN order in the 
waveform expansion given below (see Eq. (|4.ip ) after the 
second time derivative is applied. This result was already 
argued in Ref. [92j . 



III. TWO-BODY DYNAMICS WITH SPIN 
EFFECTS THROUGH 2PN ORDER 

The equations of motion in harmonic coordinates for 
the relative orbital separation a; = r n in the center of 
mass frame are USM 



1 



1 



"Newt + ~2°^1PN + Ts^SO 



1 



where 



Gm 

CtNewt = I 



-I [asiS2 + "s^ + a2PN] , (3-la) 



(3.1b) 
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OlPN 



Gm 



2 r 



n 



2f(2 - iy)v 



(3.1c) 



with m = mi + m2, v = mi m2/m^, n — x/r and t> = 
dx/dt. The 2PN acceleration given, e.g., in Ref. [5l| will 
not be needed for our calculation. The spin-orbit terms 
are (Ell 



aso = {6 [(n x v) ■ {2S + n 



(3.1d) 



-[vx {7S + 35-S)]+3r[n x (SS' + iJS)]} , 

where we denote with S = {mi — 1112) /m and 

S=Si + S2, (3.2a) 
5^2 Si 



S = 

m2 mi 

The spini-spin2 interaction terms are [sH ] 
3G 



(3.2b) 



aSiS2 



[(Si •S'2)-5(n.S'i)(n.S'2)]n 



(n • Si)S2 + (n • ^2)^1 



(3.3a) 



As originally computed in Ref. |59[ (see Eq. (|2.16p above), 
an additional term due to the influence of the spin- 
induced mass quadrupolc moment on the motion arises 
at 2PN order: 



3G 



-^Sl + q K2S2 



— n 



— in-Si)Si + qK2in-S2)S2 

q 

^{n- Sif + bqK2{n- S2f |. (3.3b) 



Here, q = mi/m2 is the mass ratio and wc recall that 
the parameters ka characterize the mass quadrupolc mo- 
ments of the bodies. 

We find that the quadratic spin contribution to the 
acceleration can be re- written in a simpler way by intro- 
ducing the spin variables 



m / Ki 



mi \K2 

m / K2 
m2 \Ki 



1/4 



1/4 



(l + Vl-KlK2)'/'Sl 
(1 - ^/l-KiK2f'^S2 , 



(3.4) 



and Sq , which is obtained by exchanging the labels 1 and 
2 in the above equation ^ Those variables generalize the 



^ In the formal limit where the induced quadrupole of at least one 
body vanishes, so that e.g. K2 — > 0, we may define the effective 
spins as: S+ = -^V2Si, S,7 = ^^Si + -^V2S2. 



quantity Sq of Ref. [60| in the case where the two bodies 
are not black holes. In terms of these spin variables the 
spin-spin part of the acceleration reads: 



3G 



[n(5o+.5„-) + (n.5„+)5o 



in ■ S^) S+ - 5n (n . S+){n ■ S^)] ■ (3.5) 



The spi n precession equations through 2PN order 
are M,m 



dS_ 

dt 



'dt 



^l^[-4{vS)-2S{vi:)]n 

+ [3{n ■ S) + S {n ■ 'E)]v + r[2S + SS] 
Gm 



(3.6a) 



-2S {v- S)- 2(1 - 2v){v ■ S)]n 
+ [Sin-S) + {l-i^)in-i:)]v 
+ f[SS+{l-2L')i:] 



(3.6b) 



It is often convenient to use a different set of spin 
variables S^^ whose magnitude remains constant and 
which obey precession equations of the form dS'^/dt = 
Ha X S'^. The relationship between the spin variables 
appearing in the equations of motion above and the con- 
stant spin variables is [4^ : 



Sr = S 



Gmv 



[2S' + (5S] 



^[{vS) + 5{v^)]v, 



Cm 

— [<5S + (l-2i.)S] 



^[6{vS) + {l-?>v){vi:)]v. 



(3.7a) 



(3.7b) 



IV. WAVEFORMS WITH SPIN EFFECTS AT 
2PN ORDER 

A. General formalism 

The gravitational radiation from the two-body system 
is calculated from symmetric trace-free radiative multi- 
pole moments 1^ and using the general formula from 
Ref. [ii truncated at 2PN order 



, TT 



2G 



1 



i?cH 3c -^"^ + 12c2^"^='^ 



gQ^3 abcde 360c^ abcdef 



cka 



1.(2) 

3c 



1 
2? 



15c3 



(4.1) 
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where N is the unit vector pointing from the center of 
mass of the source to the observer's location and R is the 
distance between the source and the observer. Here, the 
superscript {n) signifies the n**^ time derivative and the 
transverse-traceless projection operator is 



TT 



(4.2) 



2PN 
L 



1 tSO , 1 tIPN 



1 



1.5PNSO 



(4.4a) 



(4.4b) 



with Vij = 5ij - N,Nj. 



The gravitational radiation (|4.ip can be re-written in 
a post-Newtonian expansion as 



L TT 



1 



I, Newt 
"ij TT 



1 



, IPN 
2 "-ij" TT 



1 



,1.5PNSO 
'y TT 



A n-ij TT 



1 



^2PNSO 



4 'Hj TT 



^ i,2PNSS 
:a % TT 



(4.3) 



The terms in the source muhipole moments that arc a 
priori needed to compute the spin-orbit waveform exactly 
at 2PN order are identified by considering their schematic 
structure, 



T-Ncwt I 



1 



IPN 



1 



rSO 



together with the scalings of Eqs. (|4.1I) and (|3.1ap . 
Specifically, the following pieces are required: {Iab,^*'Y^^ 
using the 1.5PN motion and (/fb°)^^^ with a^''^\ from 
(Jab)'"^^ with the IPN motion and the spin evolution, 
(ji^5PNSO)(2) ^i^i^ „Ncwt^ g^^^ {J^r')^"^^ with the 1.5PN 
accurate motion, and (Ja^cd)^"^^ ^it^ a^™*-. For the SS 
part, we need {la™^)^'^^ with a^^, as the time derivative 
of /ff does not contribute at 2PN order. When we write 
the waveform in terms of the constant magnitude spin 
variables, there is an additional contribution to the 2PN 
spin piece of the waveform coming from jf^ with a^^"^^ 
and the IPN conversion factor in E'^. The relevant spin 
contributions to the multipole moments are (soj 



'1 r-\ 2 

"i 16 \ . 3 



7 7 



7 



27 109 

14 ^ 



11 47 \ 11 

4 31 \ 2Q 

7 14 7 ^ ' 14 ^ ' 



28 28 
1 Gm 2 



27 w 
(v 

14 ^ 


S) 


r n 




Gm 1 


' 29 


143 
-28" 


— +( 


-28 + 



U r 7 



(4.5a) 



Tspin 



r2 <! -^^n('nJ(i; x S)''^> " ^ - \\v)n'^'n^{v x S)'=> 



+ 3 5n<''z;^'(n X S')''^ + 3 (1 - 3i^) n<''v^'(n x Y.f^ I, 



J|Pj;' = {,5n<VJn'=5'> + (1 - 3i^)n<V^n'=S'>} . 



(4.5b) 
(4.5c) 



The nonspinning contributions to the multipole mo- Jij = —niiyr'^5eab{in-''^n°'v^- (4.6c) 

mcnts that we employed in our calculation are g Spin-orbit effects 

I-ij = mvr^ n^^n^\ (4.6a) 

Using the multipole moments of Eqs. (|4.5p and (|4.6p in 
Eq. (|4.ip and substituting the equations of motion (|3.ip 
^ijfc = —mvr^ Sn^^n-'n''^ , (4.6b) and (j3.3bp . we find the following 2PN spin-orbit piece: 
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7 2PNSO 
"-ij TT 



-(3 - Wi^) f2 (n X £c) • iV + 30(1 - Ai^){n • iV) f (n x v) • Sc 



-(7 - 29v) f (t; X Sc) • TV - 6(1 - 4:iy){v ■ N){n x • Sc - -(3 - 13i/) (n x T,^) ■ N 



2Gm, 



35 



7 



(1 - 5iy)in xT,c)-N + S[—r^{nxSc)-N--v^nxSc)-N+ 60(n • TV) r{n x v) ■ S, 



- 12{v ■ N){n xv)-Sc-l3r{vx Sc)-N 
+ 7i°(n X TV)** 



+ n"(n X Sc)''S 



35(n • TV) _ ^^4^^ • TV) f - 7(n • TV) v'^ 



^(3 - 13z^) f2 (n . Sc) - i(3 - 131.) (n • £e) + y (1 - 3z/) (n • TV)(Ar • Se) 



3 2Cm 
-5(1 - 3;.) ?• (t; • TV)(TV • Sc) - -(1 - 3;^) (n • TV)(TV • Sc) ^(1 - 3u){n ■ TV)(TV • Sc) 



+ (1 - 5zy)(n • - (3 + lliy) r (v ■ i:^) + S 

3r 



/ 4:Gm . ^ . 35 , ^ , 7 , , „ , 

(n • 5c) + — r2 (n ■ S,) - - (n • Sc) 

\ r 2 2 



{n-N){N-S, 



2Cm 3 

(n • TV)(TV ■S^-- in ■ N){N ■Sc)-5f{v- TV)(TV ■ Sc) + r (v ■ S^) 

r 2 



+ n°(n X Sc)** 
4Gm 



5(3 - liv){n ■ TV) f2 - (3 - 13i/)(n • TV) - 2(3 - lA:v){v ■ TV) f 



3r- 

(n X TV)"^!:!! 



(1 - bv){n ■ TV) 



n''{n X vf r [2(1 - 4i/)(TV • S^) + 65 (TV • S^)\ 



^(1 + 7i.) f2 + ^(1 _ 3j.)(n • TV)2 f2 _ 5(1 _ 3j,)(„ ■ N){v ■ N)r + ^(1 - 3i.)(t; • TV)^ 



+ ^(11 - 25i/)«2 - ^(1 - 3z.)(n • Ar)2 _ ^(u + 2j.) - ^(1 - 3z.)(n • Ar)^ 



(n X TV)"^^^ (5 



+ y (n • TV)2 - 5(n • Ar)(t; . Af) f + ^(t; • TV)^ + iv^ 



l{n-Nfv^~—in-Ny 
4 r 



(n X t;)''E;; (1-4:.) 



2(v • TV) - 2(n • TV)f 



36(-l + 4i/)(n • TV)(n x v) • Sc - 4(2 - 9v) f (n x Sc) • TV + -(13 - 55^)(?; x Sc) • TV 
-S ( -72(n • TV)(n x v) • Sc - 20 f (n x Sc) • TV + ^{v x Sc) • TV 
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-6(n • TV)r H (t> • TV) 

3 



+ n''(v X ScTS 



2(-7 + 29^) f (n • TV) + -(10 - 43j.)(t; • TV) 
3 

-2(5 - 22i/) f (n • TV) + - (1 - 6i/) {v ■ TV) 
3 



+ (n X v)''S^J 
+ n°(« X Sc)' 
+ i;"(n X Sc)*" 
+ i;"(n X Sc)***? 
+ z;"(n X TV)' 

+ ^(1 - 3i^){v ■ N){N • Sc) + 5 ( ^(-u • N){N ■ 5c) - 11 f (n • 5c) - 5 f (n • N)iN ■ 5c) 



-26 r(n- TV) + -(7; -TV) 



4- w°(n X 7;)' 



26f(« ■ N) + 12{v ■ TV) 



v'^iv X 5c)'5y(n.TV) 

i;"(i; X Sc)' -(16 - 67z/)(n • TV) 
3 

2(-l + 4i.)(Ar.Sc)-y(5(TV-5c) 



-(3 - 23^) r {n ■ Sc) - 5(1 - r {n ■ TV)(TV • Sc) + -(1 + 8i^){v ■ Sc) 
10 



S^{v X TV) '(5 



5 5 in 
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X AT)'' 



+ v''{v X TV)'' 



-^(1 + iy)r-l{l- f (n • Nf + ^(1 - 3iy){v ■ N){n ■ N) 



^ - iiA (n X Sc) • AT + (n x S^) ■ N 

O / o 



(Sc X Ar)°?i'' 



-(1 + 7i^)r^ + —(1 - 3i^)r2(„ . jy)2 + 5(_i + 3t/)f(n • Ar)(v • AT) + -(1 - 3i^){v ■ N)^ 



1 

T2 



^(11 - 25i/)f2 + + 3i^)(n • Ar)2w2 + ^(-17 + lOi/) + ^(-1 + 3i^)(n • Ar)2 



Cm , 
r 



(Sc X Ar)^7,'',5 
2Gm Gm 



+ • - ■ N){v ■ N) + ^{v ■ Nf + ^v^ - ^t>2(n • AT) 



5 -2 , 15 -2/ 



r r 



(Se X N^V' 

{Sc X Af)^>''5 



(n • Ar)2 

_ ^(1 + + -(-1 + 3i/)r(n • TV)^ + —(1 - 3u){n ■ N){v ■ N) 
6 2 3 

' ^ ^ 10 
-f--r{n-N)' + -{n-N){vN) 



+ {vx Ar)"n'' 



(-3 + 2'iiy)r{n ■ Sc) + 5(-l + 3i/)f(n • Ar)(Sc • AT) + -(5 + 7iy){v ■ Sc 



y (1 - 3zy)(t; • Ar)(Sc • AT) + 5 ('-llr(n • S^) - 5f(n • Ar)(5c • AT) + ^{v . Sc) + y (-y • Ar)(5c 



AT) 



(4.7) 



r 



These contributions add linearly to the other known 
terms in the waveform. Note that in Eq. (|4.7p we have al- 
ready anticipated the transverse-traceless projection and 



simplified the expression (|4.7p using that SJ^j, 
and the interchange identity (5l| 



A-iB X AT)^ = ^.(^ ^ 

for any vectors A and B. 



C. Spin-spin effects 



(4.8) 



We notice that the spin-orbit contributions at 2PN order 
arc zero for an equal-mass, equal-spin black-hole binary. 
This is a consequence of the current multipoles (|4.5p be- 
ing zero for this highly symmetric binary configuration. 

The general results (|4.7p and (|4.9p are available as a 
Mathematica notebook upon request to be used to com- 
pute the gravitational polarizations and spherical har- 
monic modes for precessing binaries for any choice of the 
source frame and the polarization triad [1, HI, [5l|, ItiI - 
[rsj . Below, we shall derive the polarizations and spin- 
weigthed spherical-harmonic modes for the case of non- 
precessing compact binaries on circular orbits. 



Spin-spin terms in the waveform at 2PN order are en- 
tirely due to the equations of motion; they arise when 
substituting a^^ in the time derivatives of I^™^- The 
second time derivative of the contribution J^j, given in 
Eq. ([2719]) is at least of 3PN order (due to the fact that 
spins arc constant at leading approximation) and there- 
fore vanishes for our calculation. We derive 



L2PNSS „ ^ 

'Hj TT - ^3fl^4 ''ijab 



5in-S+)in-So)-iS+-Sa) 



n" 5*0+ (n -So)- n" [n ■ 



D. Reduction to quasi-circular orbits 



(4.9) 



We now specialize Eqs. (|4.7p and (|4.9p to the case of or- 
bits which have a constant separation r in the absence of 
radiation reaction and for which the precession timescale 
is much longer than an orbital period. The details of 
the derivation of the modified Kepler law relating the 
orbit-averaged orbital angular frequency w and the orbit- 
averaged orbital separation are discussed in Ref. [9^ . 
The instantaneous accelerations p.ip and (j3.3bp are pro- 
jected onto a triad consisting of the following unit vec- 
tors: n = x/r, the vector £ = Ln/I^nI orthogonal to 
the instantaneous orbital plane, where = 'tij/ x x v 
denotes the Newtonian orbital angular momentum, and 
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A = {. X n. The orbital separation r and angular 
frequency to are decomposed into their orbit averaged 
piece, indicated by an overbar, and remaining fluctuat- 
ing pieces, r = f + 5r and uj = uj + Suj. Projecting the 
equations of motion along A yields the equality 2w f + ti; r 
or, equivalently 



3G 



2muj r^d^ dt 



{n-S+){n-S^). (4.10) 



At the 2PN order, r and w can be replaced by the con- 
stants r and TD respectively on the right-hand side. The 
expression for ujr^ follows from (i) dropping the time 
derivatives in the above equation, and (ii) adding an in- 
tegration constant determined by averaging w over an 
orbit. Inserting the result in the projection along n of 
the equations of motion. 



{n ■ a) 



(4.11) 



and linearizing in 5r wc find an explicit solution to the 
differential equation given by 

d5r 



2m^rc' 



2 r-{n-a) 

UJ = 



j[in-S+)iX-S-) + iX-S+)in-S-)] 



(4.12a) 



Gm r / , ^ Gm 
rc^ 

Gm\i5{£-Sc) + 3S (£• Se 



V rc 



mrc^ 



1 



-5in-S+)in-So) 



{S+ ■S^)+2{£-S+){£-S^) 

(4.12b) 



Inverting Eq. (j4.12bp to write r as a function of uj in Eq. 
(|4.7p and inserting there the expression (|4.12ap of f. we 
obtain the following spin-orbit terms in the waveform: 



h: 



2PNSO 
ij TT 



G vmuj 



'Pliab{ n^n" [4(1 - lv){t ■ Sc)(A • N) - (13 - h9>v){n x Sc) • iV - 215 (n x 5c) • N\ 



-A" A'' 
-A" 



Sc) ■ N - 262{n ■ N){£ ■ 



-FA" £'' 



3Rc^ 

4(7 - 24iy){£ ■ £'=)(A • N) + 4(1 - 6i/)(n x T,") ■ N + 5^4(n x S") ■ N + 52{£ ■ S"=)(A • N) 
4(13 - 55i/)(A X Sc) • AT + 2(-63 + 239i')(n • N){£ ■ Sc) 

f E;f^''12(l-4zy)(A-7V) 
12(-1 + Aiy){N ■ Sc) + 8(1 - 6i^)(A • Sc)(A • AT) + 4(-16 + 67iy){n ■ Sc)(n • N) 



+71" £^ 



28(7V • Sc) + 8(A • S^){X ■ N) - 128(n • S^){n ■ N) 



2(-13 + 59i/)(A • Sc)(n • AT) + 4(-10 + 43i/)(n • Sc)(A • AT) 



42(A • S^){n ■ N) - 72(n • S^){\ ■ N) 



e" 286 (A • AT) 



-hn°(n X AT)'' 
+A°(n X Nf 
+S°(n X AT)** 
+S^{n X NfS 



-(1 + iy){n ■ Sc) - 21(1 - 3i/)(n • Ar)(Ar • Sc) + (5(^3(n • Sc) - 21(n • Ar)(Ar • S, 
2(7 + 23i^)(A • Sc) + 40(1 - 3iy){N ■ Sc)(A • N) + 6(^A0{N ■ S'c)(A • TV) - 2(A • S^) 
- (21 + 17i^) + 20(1 - 3iy){\ ■ Nf + 21(-1 + 3iy){n ■ N)^ 



- 9 + 20{X ■ Nf - 21{n ■ NY' 



5c" (Ax Ar)M0(5(A- Ar)(n- AT) 
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(A X TV)' 



-80i^{n ■ Sc) + 20(1 - 3i^)(n • N){N ■ Sc) + 5 40(n • Sc) + 2Q{N ■ Sc){n ■ N) 



^E^ (A X TV)'' 40(1 - 3i^)(A • TV)(n • TV) 



(4.13) 



r 



Here, we have used that 



(4.14) 



and similarly for Sc. 

Finally, we derive the 2PN spin-spin terms for circular 
orbits. They read 



7 2PNSS _ 2Gl^W^ -prp J 

" mRc^ ' 'J°'' 
2 



40, 



+ A" A' 



-ie-S+){e-So) + -in-S+){n-S^ 



-^{n-S+){n-S^ 



211" X" 



10. 

y 

in-S+)iX-Sa) + in-S,){X-S+ 



3{n ■ S+) n^" So"' - 3(n • 5(7) n^" 



(a 



(4.15) 



E. Polarizations for non-precessing, spinning 
compact bodies 

The two polarization states and hx are obtained by 
choosing a coordinate system and taking linear combina- 
tions of the components of hj^^ . Using an orthonormal 
triad consisting of TV and two polarization vectors P and 
Q, the polarizations arc 



1 



2 yP'P^ ^Q'Q^)K^ 



TT 



TT 



(4.16) 
(4.17) 



Although different choices of P and Q give different po- 
larizations, the particular linear combination of /i+ and 
hx corresponding to the physical strain measured in a 
detector is independent of the convention used. For non- 
spinning binaries, one usually chooses a coordinate sys- 
tem such that the orbital plane lies in the x-y plane, and 
the direction of gravitational-wave propagation, TV, is in 
the x-z plane. 

When the spins of the bodies are aligned or anti- 
aligned with the orbital angular momentum, the system's 
evolution is qualitatively similar to the case of nonspin- 
ning bodies. This case is characterized by the absence 
of precession of the spins and orbital angular momentum 
and thus the orbital plane remains fixed in space. How- 
ever, the effect of the spins gives a contribution to the 
phase and a correction to the amplitude of the waveform, 
which we explicitly provide in this subsection. We use the 
conventions that the z-axis coincides with £ and the vec- 
tors n, A and TV have the following (x, y, z)-components: 



£=(0,0,1), TV = (sine*, 0, cos 61) , (4.18a) 

n = (sin$, - cos<I>, 0), A = (cos <&, sin<i>, 0),(4.18b) 

where <i> is the orbital phase defined such that at the 
initial time, n points in the x— direction. We use the 
following polarization vectors: 



P = Nx£, Q = N X P. 



(4.19) 



The vector P is the ascending node where the orbital 
separation vector crosses the plane of the sky from below. 
With these conventions, Eqs. (|4.16p with Eqs. (|4.13p . 
specialized to the case that the only non-vanishing spin 
components are (S*^ • £) and {S"^ ■ £), become: 



^2PNspin ^ _^:^cos$ sinfl {3(5 (£• Sc)(-33 + cos^ 0) + [(-93 + 167J.) + 9(1-31.) cos^e] (£-£c)} 

9 

■ cos(3$) sin 6* {5 (5 - cos^ e){£ ■ Sc) + 3(1 - Si^) sin^ 0{£ ■ Sc)} 



4i?c4 

cos(2ci>) (1 + cos^ e) {£ ■ S+)i£ -S^), (4.20) 

,2 



, 2PN spin ^ sin$ sin(20) {66 {£ ■ Sc) (-33 + cos^ 9) + [(-171 + 289i.) + 3(1 - 3iy) cos(20)] {£ ■ S^)} 

■ sin(3$) sin(26') {d {£ ■ Sc) {7 - 3 cos^ 9) + 3(1 - Sv) sin^ e{£ • Sc) } 



48i?c4 
9G^j.TOw^ 
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' mRc^ 



sin(2$) cose{e- S^){£- Sq). 



(4.21) 



F. Gravitational modes for non-precessing, 
spinning compact bodies 

The gravitational wave modes are obtained by expand- 
ing the complex polarization 

h^h+~ihy,, (4.22) 

into spin- weighted s = — 2 spherical harmonics as 



+00 e 

h{e, 0) = ^ ^ him -2Y 

1=1 m=-e 



(4.23) 



where 



/ 2/ 4- 1 

sY'^{e,<p) = (-l)y_^di„(0)e™^ (4, 



24) 



with 



mm{t+m,f.-s) ^ 



E 



/c— max(0,m — s) 



fc! 



V^(£ 4- m)\{l - m)!(£ + s)!(^ ~\ 

^ {k-m + sy.ie + TO - fc)!(^ -k~sy. 

X (cos(0/2)) (sm(0/2)) . 

The modes ft,^™ can be extracted by computing 



h,r,^ - / dn h{e, 0) _2Y 



(4.25) 



(4.26) 



where the integration is over the solid angle / dfi = 

Jq sin 9d9 J^^ d(f> and using the orthogonality property 

Jdn _^y^'"(6',(/))_^r^'™'*(6',(/)) = (5"' (J™™', where 5"' 
is the Kronecker symbol and the star denotes complex 
conjugation. Using Eqs. and in Eq. 

we find the following non- vanishing modes: 



hii = 

^33 = 



2G^mvuj^ /167r 



(4.27) 



-^<5(£.Se) + ^(-79 + 139i.)(£.Se), 



(^■5o+)(£-5o-) 

Gto2 



<^ (£ • Sc) 



(4.28a) 
(4.28b) 




24V14 



(4.28c) 



3^^(^-^c)-^^^(l-3.)(£.E.). 



hii = 7:^(5 (^•5e) + -^(l-3i^)(£-Sc) 



168^2 



I68V2 



= i^J(£.5c) + ^(l-3^)(^.Sc) 
8V14 8V14 



(4.28d) 
(4.28e) 
(4.28f) 



We have explicitly checked that in the test-mass limit 
1/ ^ 0, Eqs. (p:^ reduce to the 2PN 0{q) and 0{q^) 
terms given in Eqs. (21) of Refs. [2^,[7^, after accounting 
for the factor of (— i)™ due to the different conventions 
for the phase origin, as explained in Rcf . • 

It is interesting to note from Eq. (|4.28b[) that in the 
non-precessing case, the dominant /122 mode contains 
only terms that are quadratic in the spin at 2PN order. 
By contrast, for precessing binaries, the 2PN spin-orbit 
terms will give a non-vanishing contribution to the 22- 
modc. 



V. CONCLUSIONS 

We have extended the knowledge of the spin terms in 
the gravitational-wave strain tensor to 2PN accuracy for 
precessing binaries. Our result includes the spin-orbit 
as well as the spini-spin2 and spinj, spin| effects. The 
quadratic-in-spin terms are entirely due to the equations 
of motion, whereas the 2PN spin-orbit terms come from 
both the corrections to the orbital dynamics and the ra- 
diation field. 

For a given choice of an orthonormal polarization 
triad and a source frame, the gravitational-wave polar- 
izations can be obtained by projecting our result for 
the gravitational-wave strain tensor given in Sees. IIVBI 
and IIV CI orthogonal to the propagation direction. For 
precessing binaries, there is no preferred unique choice of 
the source frame [1, [s^ [Sll, [7ll - [75j , but in the case that 
the spins are coUinear with the orbital angular momen- 
tum, the procedure to obtain the polarizations can be 
carried out in a similar fashion as for nonspinning bina- 
ries. For the non-precessing case and circular orbits, we 
provided rcady-to-use expressions for the gravitational 
polarizations in Sec. IIV El which could be directly em- 
ployed in time-domain post-Newtonian, phenomenolog- 
ical and effective-one-body-based template models |19l - 
[2l,[5l|,[9|. 

In view of the current interest in interfacing analyt- 
ical and numerical relativity, we also provided the de- 
composition of the waveform into spin- weighted spherical 
harmonic modes for non-precessing binaries and quasi- 
circular orbits. We verified that the test-particle limit 
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of our result reduces to the expressions obtained from 
black- hole perturbation theory [tI, [tI] . We noted that 
for spins collinear with the orbital angular momentum, 
the dominant /122 mode of the waveform contains only 
quadratic-in-spin effects since the spin-orbit contribu- 
tions vanish in this case, although they are non-zero for 
generic, processing configurations. 
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Appendix A: Useful identities 

According to the way the waveform is computed, the 
result may take various forms, which arc not immediately 
seen to be equivalent. Their difference vanishes due to 



some dimensional identities valid in 3 dimensions. They 
all amount to expressing the fact that a tensor with four 
antisymmetrized indices must vanish. We shall present 
here two of such identities, which turned out to be par- 
ticularly useful for our checks, together with Eqs. (5.2) 

of Rcf. m. 

Let Ua = U\, for A e {1, 2, 3}, be three vectors of K^. 
The first identity tells us that for any vector [/, we must 
have 

(C/i X U2)^nUp{Ui ■ U) - Ui>{U3 ■ U)] (Al) 

= ui'[{u X UiyHu2 ■ c/3) -{ux U2y\ui ■ c/3)] 

+ C/^'[(J7 X U2y\Ui ■ Ui) -{Ux Uiy\U2 ■ C/4)] . 

To show this, we compute e^ab'^™^''^mpqUiU2U^Ul in two 
different manners: (i) we group the first two epsilons, 
which are next expanded in terms of the identity tensor 
S^j using the standard formula Ciafce™^*'' = •^^■^"[i^^a^\]i 
(ii) we group the last two epsilons and apply the con- 
tracted version of the previous equation: e'^^'^empq = 
2(5-'jp(5'^^j . One of the remaining free indices, say k, is fi- 
nally contracted with Uk- 
The second identity reads: 



5'^[Ululul - Ul{U2 ■ U^f - UliUz ■ Uif - UliUi ■ U2f + 2{Ui ■ U2){U2 ■ U^yu^ ■ Ui)] 
+ 2ui'u(>[U^(U3 ■ Ui) - [Ui ■ U2yU2 ■ c/3)] + 2u['U'^[Ul{Ui ■ U2) - (U2 ■ U^yu^ ■ Ui)] 

+ 2c/^V^')[c/f (c/2 • c/3) - (c/i • c/2)(c/i • c/3)] + uiui[{U2 ■ u^f - uiul] + uiuiiiu, ■ Usf - ului] 

+ UlUi [(C/i • C/2)2 - UlUi] = . (A2) 



It is proved by contracting the equality C/j''[/|C/|J']-' = hand side of Eq. (|A2p is identically zero, the non-trivial 
with UiaU2bU3c and expanding. As the trace of the left- content of this identity consists of its STF part. 
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